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We study three practical implementations of the Overlap-Dirac operator B = 
^[1 + 75e(H w )] in four dimensions. Two implementations are based on different 
representations of e(H„,) as a sum over poles. One of them is a polar decomposition 
and the other is an optimal fit to a ratio of polynomials. The third one is obtained by 
representing e(H m ) using Gegenbauer polynomials and is referred to as the fractional 
inverse method. After presenting some spectral properties of the Hermitian operator 
, we study its spectrum in a smooth SU(2) instanton background with the 
aim of comparing the three implementations of B D . We also present some results 
in SU(2) gauge field backgrounds generated at j3 = 2.5 on an 8 4 lattice. Chiral 
properties have been numerically verified. 
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1 Introduction 

The overlap formalism provides a way of realizing exact chiral symmetry on 
the lattice. For vector gauge theories, there is now promise that this formalism 
can be made practical. The Overlap-Dirac operator derived from the overlap 
formalism [1] of chiral fermions is [2] 
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is the Hermitian Wilson-Dirac operator. We will refer to m as the overlap 
mass. This is a parameter that has to be in the range (0, 2) for B G to describe 
a single massless Dirac fermion. In principle, any value of m in this region 
will yield the same continuum theory. But at finite lattice spacing the cutoff 
effects can be quite different, for different choices of m. In particular one needs 
m > mi(g 2 ), for some m^g 2 ) going to zero in the continuum limit, in order 
for the overlap fermions to "feel the topology of background gauge fields" [4] . 

It is necessary to efficiently deal with the action of e(H w ) on a vector for prac- 
tical implementations of this operator in four dimensional theories. Recently, 
one representation of e(x) as a sum of poles of x 2 was used to demonstrate 
that a practical implementation of the Overlap-Dirac operator is possible in 
three dimensions [5]. One can also represent e(x) as an optimal rational func- 
tion [6]. Another approach is to use Gegenbauer polynomials to represent ^= 
resulting in an iterative procedure to compute the action of e(R w ) on a vec- 
tor [7]. In this paper we compare these three methods from a practical point 
of view. We start by presenting some spectral properties of H D = 7sp and 
~p and then use a smooth SU(2) instanton background as an example to il- 
lustrate these spectral properties. The example will also enable us to compare 
the three practical implementations of e(H w ). We also present some results in 
SU(2) gauge field backgrounds at (5 — 2.5 on an 8 4 lattice. 



2 Spectral properties of o 



In this section, we build on some results of Neuberger [2,3] by considering the 
eigenvalue problem of the Hermitian Overlap-Dirac operator, 

H = 75 B = i[ 75 + e(H w )] . (5) 

We do not use the explicit expression for H w in this section and we only use 
the fact that e 2 (H^) = 1. The exact chiral symmetry built into the overlap 
formalism has been reemphasized in the context of the fermionic action with 
p in Ref. [8]. We therefore expect the spectrum of H to essentially repro- 
duce all the important continuum properties. The spectral properties of the 
continuum operator H = 75D are 



2 



• Exact zero eigenvalues of H are associated with topology. The associated 
eigenvectors are also eigenvectors of 75. The zero eigenvalues are not paired 
in any sense. We can have n + zero eigenvalues with positive chirality and 
n_ zero eigenvalues with negative chirality. The difference (n + — nJ) is the 
topology of the background gauge field. 

• Non-zero eigenvalues of H come in pairs that are equal in magnitude and 
opposite in sign. The associated eigenvectors are not eigenvectors of 75, but 
rather 75 has zero expectation value in the eigenvectors, = 0- 

On a finite four dimensional lattice, the matrix H D is a finite even dimensional 
matrix and det H Q is the overlap formula in a fixed gauge background. The 
overlap exactly vanishes if the topology is non-trivial; therefore, H D should 
have exact zero modes in a gauge field background that is topologically non- 
trivial. These zero modes should have definite chirality and the difference in 
the positive and negative modes is the topology of the gauge field. If the 
topology is odd, then H D should have an odd number of exact zero modes. 
Therefore it should also have an odd number of non-zero modes. This im- 
plies that the number of positive eigenvalues cannot be equal to the number 
of negative eigenvalues. It is then interesting to explore the properties of the 
spectrum of H D with the aim of establishing the equivalence of the two contin- 
uum properties listed above. We prove the following properties of H Q on the 
lattice: 

• The spectrum of H G is bounded by the region [—1,1]. 

• Zero eigenvalues of H G have a definite chirality and need not occur in pairs. 
Eigenvalues of H G equal to ±1 also have definite chirality. These eigenvalues 
also need not occur in pairs. 

• Non-zero eigenvalues of H D that have a magnitude less than one come in pairs 
that are equal in magnitude and opposite in sign, namely ±A. The associated 
eigenvectors are not eigenvectors of 75, but rather 75 has expectation value 
±A in the eigenvectors, ^^5^ = ±A for R ip = ±\ip. 

Since H G is an even dimensional matrix, the unpaired zero eigenvalues have to 
be matched by unpaired eigenvalues equal to ±1. This is what is expected to 
happen in a topologically non-trivial background. 

We present the proofs to the properties of H listed above. 

• Squaring H in (5) we get 



R 2 = -[2 + 7 5e (H u ,) + e (H w ) 75 



(6) 



The above equation can be rewritten using (5) as 



2R 2 = 75 H + H o75 



(7) 
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and is commonly referred to as the Ginsparg- Wilson relation. As pointed out 
in [2], 7 5 e(H^) is an unitary operator and e(H„,)7 5 is its Hermitian conjugate. 
Therefore, all the eigenvalues of [7 5 e(H w ) + e(H w )7 5 ] are bounded by [—2, 2]. 
Therefore all eigenvalues of are bounded by [0, 1] and from this the first 
property of H c follows. 

• Let us first find all eigenvectors of H G that are also eigenvectors of 75. Let ip 
be an eigenvector of H G with eigenvalue A and let 75^ = With the aid 
of (7) one concludes that A(A =F 1) = 0. Therefore we have either A = or 
A = ±1, and we have shown that chiral eigenvectors of H have eigenvalues 
equal to zero or ±1. The subspace spanned by the zero eigenvalues of H D 
can be split into two further subspaces. One with dimension of n + will have 
positive chirality and one with dimension of n_ will have negative chirality. 
The number n + need not be equal to n_ and if (n + — nJ) is not equal to 
zero then the gauge field carries a non-trivial topology. All eigenvalues of 
H D equal to +1 have positive chirality and vice- versa. 

• Let < A < 1 be an eigenvalue of H D with eigenvector ip that is normalized 
to unity. Clearly, 

^l^ZMhM; ^0 = (8) 

is a new vector that is orthonormal to ip since ip^'juip 7^ ±1. We can now 
compute H o using the Ginsparg- Wilson relation, (7), and find that 

H ^ = - A0+ 2 ^-^) > ( 9 ) 

This shows that H D acting on the subspace spanned by (ip, 0) results in 
vectors in that subspace. From B. ip = \ip and (flip = 0, it follows that 
^Hotp = 0. From the hermiticity of H D , we also have t/^Ho^ = 0. This along 
with (9) results in 

V>W = A (10) 

since we have assumed A 7^ 0. We have also assumed that A 7^ 1 in the 
beginning and therefore is indeed a non-zero vector with 

H o 0=-A* *=^f§f ■ (ID 

This establishes the last property of H . 

Comparing the properties of H c with the desired properties of H in the con- 
tinuum, we conclude that H reproduces topology correctly, as expected. The 
eigenvectors of the non-zero eigenvalues of H D do not have a zero expectation 
value of 75. Eigenvalues of H G have to be scaled by the lattice spacing to get 
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the eigenvalues of the continuum operator and then only the eigenvalues of H G 
close to zero have a continuum relevance. The last property tells us that the 
corresponding eigenvectors also have an expectation value of 75 close to zero 
as expected in the continuum. 

It is trivial to find the eigenvalues and left eigenvectors of D G since D also maps 
the subspace spanned by (?/>, <ft) onto itself. The eigenvalues of D are zero when 
A = 0, unity when A = ±1 and the complex conjugate pairs (A 2 ± — A 2 ) 
for < A < 1. This proves that the only real eigenvalues of the Overlap-Dirac 
operator are zero and unity. The massive Overlap-Dirac operator is given by 
po+rrif with —1/2 < rrif < 00 [9] Q. Clearly the operator has a zero eigenvalue 
only when rrif — 0. Occurrence of zero eigenvalues for a lattice Dirac operator 
at positive quark masses in certain gauge field backgrounds are referred to as 
"exceptional" configurations. Such "exceptional" configurations are actually 
generic for the Wilson-Dirac operator but non-existent for the Overlap-Dirac 
operator. As pointed out in Ref. [10] the only difficulty arises in defining the 
Overlap-Dirac operator if has a zero eigenvalue since e(H„,) is then not 
well defined. In such a background configuration, one can define the operator 
at a value of m slightly away from the one where H w has a zero eigenvalue 
since the continuum limit is expected to be independent of m. 

Finally, one can also trivially get the spectral flow of the massive H (m/). 
If H o (0) has a zero eigenvalue then H (m/) has an eigenvalue equal to ±m/ 
corresponding to positive and negative chirality of the eigenvector. If H o (0) has 
an eigenvalue equal to ±1, then H D (m/) has an eigenvalue equal to ±(1 +mj). 
For the paired eigenvalues ±A of H o (0) in the range < A < 1, H D (mj) has 
paired eigenvalues equal to ±^A 2 + 2X 2 rrif + mj. This again shows that zero 
crossings in the region —1/2 < rrif < 00 occur only at rrif = and these are 
due to the topological content of the gauge field. 



3 Algorithms 

The key element in the application of the Overlap-Dirac operator (1) on a vec- 
tor is how to apply e(H w ) on a vector. The meaning of e(H w ) is defined through 
the e function on the eigenvalues of H w . Hence, the diagonalization of H„, is 
one possible method and has been applied in two dimensional systems [11]. 
However, on larger systems of interest, exact diagonalization is prohibitive. 



1 The fermion mass rrif should be written as rrif = j^— with fi in the range [—1,1]. 
Positive values of ft correspond to positive quark masses and vice-versa. [9] 
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3.1 Fractional inverse method 



By representing e(H TO ) as 

e(E w ) = (El)' 112 H w (12) 

we can directly apply e(H w ) on a vector by using a Krylov space iterative 
solver for the fractional inverse based on Gegenbauer polynomials [7]. In this 
method, one parameterizes the system to be solved, 

M> = X , (13) 

by 

M = c(l + t 2 - 2tA) , (14) 

The key to solving (13) is then the generating function for Gegenbauer poly- 
nomials 

oo 

(l+t 2 -2tAr< = J2t n C2(A) . (15) 

n=0 

Using the recursion relations for the C% one can construct an iterative solver 
for in (13). We refer the reader to [7] for details. For the case of interest, 
7 = 1/2 and M = is hermitian and positive (semi-)definite. We restrict the 
bounds of A to be in [—1, 1] and obtain uniform convergence of A and hence 
M = in (15) provided \t\ < 1. It is important for the best convergence to 
match the extremal eigenvalues of e [A^ in , A^ a J to those of A e [—1,1]. 
Then, the t and c are determined and 

t = l--, (16) 



where k = |A max /A m j„| is the condition number of H w . The parameter t is 
the convergence factor for the algorithm. Typically, \ ma x & 8, but |A m j n | can 
be quite small [4]. Control of the lowest mode is important since the number 
of iterations expected to reach a given accuracy R is bounded by — \og(R)n, 
and numerical tests have shown this bound to be saturated. We remark that 
the same bound also applies to the Conjugate Gradient (CG) algorithm but is 
not expected to be saturated. The polynomials generated by the Gegenbauer 
expansion are not optimal, while the CG algorithm is optimal. Eigenvectors 
with small eigenvalue can be projected out of (12) to improve convergence. 
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3.2 Pole methods 



3.2.1 Polar decomposition 

Another approach is to consider e(z) as the limit of some function that switches 
quickly between ±1 on z changing sign. Recently, Neuberger showed [5,12] that 

e{z) = lim f N (z), (17) 

N—HX3 



with 



f (z) - + 



(1 - z) 2N 



i + z y N + (i - z y N 



(18) 



In particular, /jv(z) provides a good approximation for e(z) at finite N for 
l/N < z < N, and can be expressed quite simply, by matching poles and 
residues, as a sum over poles 

N 

M*) = *E^ (is) 

k=l z ^~ u k 

" = s4(R»' k=h, ">-5» ' (20) 



Since e(sz) = e(z) for positive s, the argument z can be rescaled into the 
region suitable for the approximation. For the application of /^(H^) on a 
vector, one can use a multi-shift CG solver [13] for the N terms in f N (R w ) [5]. 
The convergence is governed by the smallest shift; however, the overall cost is 
not much more than the inversion for the smallest shift. The representation 
of /w(H w ) in terms of poles is the key observation to making the algorithm 
practical. 

While the smallest shift, b±, in (20) governs convergence, the choice of iV is 
determined by the condition number of H w . The optimal choice for the scale 
factor s is , 1 The maximum deviation of f^isz) away from unity 

V ^ m a x A rn, i n 

occurs at z — \ ma x and z = A min , and the maximum deviation monotonically 
decreases in N. The desired accuracy 5 over the spectral range of H w fixes iV 
depending on the condition number k, namely iV = 0.25-^/Km |. Hence, the 
condition number is the real measure of convergence for small \ m in if a good 
uniform approximation is to be achieved. This requirement may not be needed 
in practice as long as the lowest modes are well represented, but will be used 
throughout this work since it is not difficult to achieve. 
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3.2.2 Optimal rational approximation 



Other approximations to e(z) over a finite interval can be suggested. A poly- 
nomial approximation like Chebyshev over the interval [—1,1] is particularly 
bad since it is only logarithmically convergent. A better choice is the ratio of 
polynomials called an optimal rational function approximation. Since e(z) is 
odd, we choose the form 



where we assume the polynomials are irreducible of degree N. We choose the 
highest order coefficient in Qn(z) to be one. The normalization is contained 
in Pn(z) and is defined to be cq. As in the previous section, since e(sz) = e(z) 
for positive s, the argument z can be rescaled into the region suitable for the 
approximation. 

To cast into a suitable form, we change variables to y = z 2 and consider 
h(y) = (l/z)e(z) over the interval [0, z^^]. One can show [14] that there exist 
best approximations of r(y) = PN(y)/Qj^(y) for a continuous function in a 
uniform norm over the interval [y m in, 1] (with no loss of generality, we restrict 
the upper bound). The deviation of r(y) from the intended function 1/^/y can 
be shown [14] to have 2 + 2N extrema with two always the endpoints. An 
algorithm, sometimes called the Remez algorithm [6,14,15], exists to find the 
best approximation r(y). The version used here requires the extrema of the 
relative error to be constant. For our application, this implies the worst error 
in <7jv(z) is for the smallest z. Other weighting functions can be chosen. 

The numerical evaluation of (21) is most stable in factored form 



where we find the roots pk and qk to be real and positive. We can cast this 
expression as a sum over poles as in Section 3.2.1: 



We choose the value of N so that the deviations of |<?jv(z)| over the interval 
[z m i n , 1] are some acceptably small value 5. We show examples in Fig. 3 of 
N = 6,8 and 10 for the fit interval [0.005, 1], which admits a condition number 
of 200. We see the method rapidly converges in N. In our tests, we have chosen 
N = 8 where we see less than 0.5% of an error. However, we observe that since 




(21) 




(22) 




(23) 
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we are approximating a flat function, we can in practice extend the acceptable 
upper bound beyond 1. For z > 1, gN{ z ) monotonically increases. For N = 8, 
we can choose z max = 2.8 and still maintain a 1% error. We can use this trick 
to effectively increase the useful allowable condition number to over 500. In 
Table 1, we show the roots and normalizations for the polynomials shown in 
Fig. 3. 



4 Results 



We will consider a gauge field background with a single SU(2) instanton to 
compare the performance of the three algorithms described in the previous 
section. The SU(2) link elements on the lattice obtained from an exact evalu- 
ation of the path ordered integral of the continuum instanton take the following 
form [16]. In the singular gauge, 



U^{n) = exp 



a(^(n;0)-^{n;p)) 



(24) 



*V(^;p) 



: tan 



^/p 2 + E^(n,-c u y P 2 + E„K - c,) 2 + (?V - c^) 



(25) 



bi(n) = (-n 4 + c 4 , n 3 - c 3 , -n 2 + c 2 ), 
M 71 ) = ( _n 3 + °3, -ri4 + c 4 , Til - Ci), 
b 3 (n) = (n 2 - c 2 , -rii + ci, -n 4 + c 4 ), 

h(ri) = (ni - ci,n 2 - c 2 ,n 3 - c 3 ) (26) 

Here c denotes the center of the instanton, and p is the size of the instanton 
measured in lattice units. We work on an 8 4 periodic lattice and set p = 1.5 
and c M = 4.5 with the lattice sites in each direction numbered from 1 to 8. 

We first present the results. We used the explicitly real representation of H„, [3]. 
We computed a few low lying eigenvalues of H c using the Ritz algorithm [17] to 
compute the lowest eigenvalues of H 2 ,. From our previous work on the spectral 
flow of the hermitian Wilson-Dirac operator, we know that a level of H„, crosses 
zero between m = 0.5 and m = 0.6 [18]. This flow is shown as a solid line in 
Fig. 1. The low lying eigenvalues of H G as a function of m are also shown in 
the same figure as octagons. For < m < 0.5, fourteen low lying eigenvalues 
are shown and for 0.6 < m < 0.9 fifteen low lying eigenvalues are shown. 
The non-zero eigenvalues come in opposite pairs with their chirality being 
equal to their eigenvalue. This is in accordance with the last property of H . 
The seven positive and the seven negative eigenvalues are nearly degenerate. 
For 0.6 < m < 0.9, there is a single zero eigenvalue for H . This is due 
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to the instanton background. The abrupt appearance of a zero eigenvalue 
with chirality +1 is accompanied with the appearance of an eigenvalue of —1 
with chirality —1. It is difficult to confirm this numerically because of the 
observed high density of eigenvalues near ±1. The shape of the zero mode for 
0.6 < m < 0.9 is shown in Fig. 2. The z(t) plotted is equal to Y.xa^i^^^ a ) 
and the continuum mode is plotted as a solid line. The zero mode is remarkably 
stable under change of the overlap mass m. 

We have made comparisons of the average number of operations to obtain 
the result of the action of e(H w ) on a vector. The average is over all the Ritz 
iterations needed in solving for the lowest eigenvalue in Fig. 1. The main 
comparison is between the pole method and the fractional inverse method. 
The two pole methods are expected to be comparable in the average number 
of operations on a vector. Since we are using the method of pole shifts 
the number of poles only affect the memory requirement and some additional 
multiplications of vectors by scalars and additions of vectors. We show in Fig. 4 
the scaling of the number of operations versus the condition number k of 
H„, at masses below and above the crossing. We used the N = 8 polynomial in 
Tab. 1 (c.f. Eq. (23)). The fractional inverse method scales quite linearly in k 
and the operation count grows to large values while in the optimal rational 
approximation it grows very slowly. This is because the bound on the rate of 
convergence is saturated in the fractional inverse method while in CG it is 
not. As expected, no substantial difference in the number of operations 
was found between a N = 50 order /at(H„,) in Eq. (20) and the N = 8 optimal 
rational approximation. 

Although the two pole methods need similar number of operations there are 
some differences. For a fixed accuracy (maximum deviation 5 away from unity 
in the realization of e(z)) and a fixed condition number the number of poles 
needed in the optimal rational approximation is smaller than the one needed 
in the polar decomposition. For example, with 5 = 0.01, the optimal rational 
approximation needs seven poles for a condition number of 200 and twelve 
poles for a condition number of 1000. In contrast, the polar decomposition 
needs 19 poles for a condition number of 200 and 42 poles for a condition 
number of 1000. This indicates that the number of poles needed grows with 
the condition number at a slower rate in the optimal rational approximation 
than in the polar decomposition. 

A comparison of the number of poles needed for the optimal rational approxi- 
mation as a function of the accuracy for a fixed condition number is shown in 
Fig. 5. The number of poles needed depends logarithmically on the accuracy, 
5, and the slope is smaller for larger condition number. As mentioned before, 
the number of poles needed only affects the memory requirement and results 
in a few additional operations of two types: addition of two vectors and mul- 
tiplication of vector by a scalar. For a fixed memory requirement, the optimal 
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rational approximation will admit a larger condition number. 

We also present the results of the spectral flow of H Q in three pure SU(2) gauge 
field backgrounds generated at (3 = 2.5 on an 8 4 lattice to further emphasize 
the feasibility of the algorithms described here. The flows are shown in Fig. 6- 
8. Fig. 6 is the flow in a background that contains a single instanton and 
this results in an exact zero mode of H D . We verified that it is an exact zero 
mode by computing its chirality which was found to be equal to —1. Fig. 7 
is the flow in a background that has no topological object and we do not 
find any exact zero mode in this case. Fig. 8 is the flow in a background that 
contains an instanton and an anti-instanton. These will result in approximate 
zero modes and one can see evidence for this in the sharp drop in the flow 
between m = 0.6 and m = 0.7 and also by comparing to the flow in Fig. 7. 
The modes in Fig. 8 are definitely not exact zero modes and, in accordance 
with the last property of H G , have chiralities equal to their eigenvalues. Finally 
we also show the exact zero mode found in the flow in Fig. 6. We plot it for 
all m > 0.8 and find that the mode is essentially the same and it also agrees 
with the zero mode of H^, computed close to the point in m where H w has a 
zero crossing. 



5 Concluding remarks 

We have compared three implementations of the Overlap-Dirac operator in 
four dimensions. We find that the pole methods are superior in convergence 
to the fractional inverse method. The main advantage of the optimal rational 
function approximation is that a suitable error can be achieved for a smaller 
order polynomial. This saves some work in the multi-shift CG since fewer 
vectors are need. Also, memory can be saved. No difference in convergence 
properties is expected otherwise. 

For the massive Overlap-Dirac operator, one can again use the multi-shift CG 
solver [13] to solve for the inverse with several quark masses. The action of 
e(H w ) is only applied on one vector and the quark masses can be expressed as 
a shift. 

In practice it may be worth projecting out a few small eigenvectors of B. w before 
the action of e(H. w ) on a vector. This is useful if there are a few well separated 
small eigenvalues. The projection can have a dramatic effect in the fractional 
inverse method since the number of operations is linear in the condition 
number. It will also decrease the operation count in the pole methods. We 
have verified these effects on the configurations studied in this paper. 

The pairing of the non-zero eigenvalues of H D can be used to reduce the com- 
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putation of the spectrum of H D . 

We have demonstrated the feasibility of using the Overlap-Dirac operator in a 
four dimensional gauge theory by computing the spectral flow of H D in SU(2) 
gauge field backgrounds at /3 = 2.5 on an 8 4 lattice. Chiral properties have 
been shown to be properly reproduced. Exact zero modes have been shown to 
be related to topology. They are clearly differentiated from almost zero modes 
in that the exact zero modes have chirality ±1 while the almost zero modes 
have chirality equal to their eigenvalue, namely very small. 
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cj, 6) = 0.182031 


c{, 8) = 0.138827 


c 10) = 0.112209 


Pe Qe 
2.05911386e-04 2.59134855e-05 
4.41581797e-03 1.09739974e-03 
3.77980526e-02 1.40770643e-02 
2.01546314e-01 9.03340757e-02 
1.00302812e+00 4.40335697e-01 
1.18344772e+01 2.67694284e+00 


^8 Qs 
8.76771341e-05 1.41842155e-05 
1.23623252e-03 3.61450316e-04 
9.11319209e-03 3.60062816e-03 
4.26239090e-02 2.05879480e-02 
1.54018532e-01 8.27820958e-02 
5.10199622e-01 2.80210832e-01 
1.96142588e+00 9.59252529e-01 
2.06291295e+01 4.86700520e+00 


Pio Qw 
4.66940445e-05 8.66147628e-06 
4.91389639e-04 1.63943675e-04 
3.13861891e-03 1.31054157e-03 
1.36974237e-02 6.82851948e-03 
4.57939021e-02 2.57134015e-02 
1.300671276-01 7.83188602e-02 
3.43184790e-01 2.12031104e-01 
9.33072953e-01 5.59010141e-01 
3.19171360e+00 1.64058378e+00 
3.17877409e+01 7.65233478e+00 



Table 1 

The normalization and roots for the optimal rational function approximation to 
e(z) (Eq. 22) over the interval [0.005, 1] for N = 6, 8 and 10. 
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0.0 0.2 0.4 0.6 0.8 

m 

Fig. 1. The octagons describe the low lying spectrum of H Q in an smooth SU(2) 
instanton background. The zero mode found at m > 0.5 is singly degenerate and is 
associated with the instanton. The lines are the spectral flow of H w . 
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Fig. 2. The mode associated with the zero eigenvalues in Fig. 1 at m = 0.6 (square), 
m = 0.7 (diamond), m = 0.8 (octagon) and m = 0.9 (plus). The solid line is the 
continuum zero mode with p = 1.5. 
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Fig. 3. The optimal rational function approximation gjv(z) to e(z) in the interval 
[0.005, 1] for N = 6, 8 and 10. Shown is a detail of the region of largest deviation. 
Outside this region, the curves rapidly approach unity. 
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Fig. 4. Average number of operations for the optimal rational approximation 
and the fractional inverse methods to apply e(H w ) on a vector versus the Wil- 
son-Dirac operator condition number for the configuration in Fig. 1. The cross and 
fancy cross are for the optimal rational approximation below and above the crossing, 
resp. The diamond and fancy diamond are for the fractional inverse method below 
and above the crossing, resp. 
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Fig. 5. The accuracy as a function of the number of poles for a fixed condition 
number. The crosses are for condition number k = 200 and the octagons are for 
k = 1000. 
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Fig. 6. The octagons describe the low lying spectrum of H G in a pure SU(2) gauge 
field background at (3 = 2.5. The zero mode found at m > 0.8 is singly degenerate 
and is associated with an instanton in the background gauge field. The lines are the 
spectral flow of H w . 



20 




Fig. 7. Same as in Fig. 6 for a different SU(2) configuration at (5 = 2.5. 
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Fig. 9. The mode associated with the zero eigenvalues in Fig. 6 for m > 0.8 
plotted with octagon symbols and compared with the zero mode of H^, shown 
straight line segments. 
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